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SUMMARY 


This paper describes a discrete modal-control design approach that is applied to 
a single-control-surface, unswept aircraft wing subject to bending-torsion flutter. 
The modal approach is a mathematical method to decouple the equations of motion into 
isolated differential equations. In this paper, a pole-placement approach is then 
applied to determine stability gains in the discrete plane using only the two 
comp lex- con jugate flutter-mode equations. A fixed-gain Kalman filter is used to 
estimate the modal amplitudes using three measurements. Results are presented for a 
full-state estimator (36 states) and two reduced-state estimators using two different 
closed- loop pole locations. The control law is designed for a dynamic pressure that 
is 50 percent greater than the uncontrolled-f lutter dynamic pressure. With constant 
control-law gains, the closed-loop system remains stable over the dynamic-pressure 
range from flutter onset to approximately an 80-percent increase in pressure. 


INTRODUCTION 

In the 1960’s and early 1970* s, research in controlling the primary mirror of a 
large orbiting space telescope was conducted. The goal was to control the figure 
error of the primary mirror to a fraction of a wavelength with many actuators mounted 
at the rear of the mirror and with a figure-error sensor located at the center of 
curvature. A modal-control technique was developed (ref. 1) to decouple the equa- 
tions of motion representing the plant. The technique allowed each controlled mode 
to be individually compensated for in the modal domain. Additional research applying 
the modal-control approach to the mirror problem is described in references 2 and 3. 
Subsequent to this work there has been renewed interest in applying modal control to 
other areas, particularly to the area related to large, flexible spacecraft control 
(refs. 4 to 7) . Another area in which modal control could have a significant impact 
relates to flutter suppression. This paper describes a discrete modal-control design 
approach that is applied to a single-control-surface, unswept aircraft wing subject 
to bending-torsion flutter. One advantage of the modal-control approach is that the 
flutter mode, as well as all other flexible-wing modes, can be mathematically repre- 
sented by two isolated complex- conjugate equations. The decoupled nature of the 
modal- control approach allows additional insight into the control problem, and it is 
possible to use either classical control techniques or modern control theory to 
design feedback gains. In this paper, a pole-placement approach is used to calculate 
the stability feedback gains, which are analytically designed by using only two 
complex-conjugate flutter-mode equations. 

Several other approaches have been applied to the flutter-suppression problem. 
These approaches include an aerodynamic energy method (refs. 8 and 9) and optimal- 
control approaches (refs. 10 to 12). A modal-control analysis is used in refer- 
ence 12, but the design is for a continuous-control law and the performance index is 
expressed in terms of a quadratic cost function that is related to all of the modal 
amplitudes. As described above, the main thrust of this paper is to isolate and 
control only the flutter mode, influencing the higher order modes as little as possi- 
ble. The modal-control approach appears to be a good method to meet this goal. 



The first section of this paper includes a brief description of the v/ing and the 
various subsystems that comprise the open-loop plant. Detailed models of all subsys- 
tems are given in appendix A. 

The second section includes the equations for transforming the physical domain 
to the modal domain and includes the form of the various matrices. A theoretical 
solution for the modal equations of motion is presented for a discrete sampling 
period. The analysis carries through to the use of the flutter-mode equation to 
design the stability feedback gains. 

While the research was being conducted, it became obvious that wind gust and the 
external disturbance states should be estimated. The approach used in the third 
section is to augment the modal equations and solve a convolution integral to deter- 
mine the effect of these states on the modal amplitudes. The discrete form for the 
total system with 36 state equations is shown. 

The fourth section of this paper shows how the discretized fixed-gain Kalman 
filter is calculated. The symbology illustrates the matrices that are to be approxi- 
mated both for off-design analysis and for reduced-state estimation. Two methods are 
presented for calculating the reduced-state estimator. Also included are the equa- 
tions to calculate the estimator closed- loop stability and the total closed-loop- 
system stability. 

The fifth section ("Results") contains several figures showing the performance 
for two design eigenvalue locations. One eigenvalue has a low damping ratio of 0.14 
and the other a higher damping ratio of 0.62. Performance is determined by normal- 
izing the root -mean- square (rms) value of each variable about a zero mean with the 
rms value of wind gust about a zero mean. The plots shown compare a base run using 
perfect modal feedback and stability gains adjusted for each dynamic pressure with 
several runs using three output measurements with an estimator. Both full-state 
estimation and reduced-state estimation runs are presented. The "Results" section 
also contains figures showing some typical data signals for several states and step 
responses due to an external command input. 


DESCRIPTION OF MODEL 

An unswept aircraft wing (fig. 1) subject to bending-torsion flatter was used as 
a model for application of the modal-control concepts. Most of the data for this 
wing are in reference 13. The semispan for this wing model is 5.08 m (200 in.) and 
the chord varies linearly from 2.54 m (100 in.) at the root to 1.52 m (60 in.) at the 
tip. Five grid points, excluding the root section, are located along the elastic 
axis at the corresponding 40-percent-chord location. The aerodynamic center is 
located at 25 percent chord. 

The wing model has the capability of incorporating one or more control surfaces 
along either the leading or the trailing edge or along both edges. For all analyses 
in this paper, a single trailing-edge control surface is used. The equations of 
motion that describe the complete system are derived from several detailed models 
described and defined in appendix A. These detailed models include the structural 
equations of motion, unsteady aerodynamics with Jones' approximation of the Wagner 
function (refs. 14 and 15), an external-disturbance model using the Kussner function 
(ref. 16) to relate the wind gust to forces and moments on the wing, a Dryden wind 
model (ref. 17) for the vertical direction, and an in-line compensator that relates 
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the feedback control signal to the actuator input signal. When all of these models 
are combined as in appendix A, the state- space equation of motion for the complete 
system is 


X = AX + BA + DX^ 
m mm me m d 


( 1 ) 


with X^ defined as 


X = 
m 


w 


w 


X 


X. 


( 2 ) 


where is the open- loop-model state matrix. 


m 


is the open-loop-model control 


is the open- loop-model disturbance input matrix, X^^ and X 


m 


IS 

and 


input matrix, 

represent the state vector and its time derivative for the open-loop model, A^ 
the control feedback signal, X^ is the external-disturbance state vector, X^ 

X represent the wing state vector and its time derivative, X^ is the unsteady 
aerodynamic state vector, and X^ is the compensator state vector. (A list of sym- 
bols used in this paper appears after the references.) A block diagram of the open- 
loop- system model is shown in figure 2. 


MODAL ANALYSIS 


Modal analysis is an approach for transforming the system (eq. (1)) into an 
equivalent number of decoupled state-space equations. The rationale for this 
approach is that the decoupled equations allow added insight into the control prob- 
lem. Stability feedback gains can be calculated with only the two complex unstable 
flutter-mode equations by either a classical control approach or modern control 
theory. Since the equations are decoupled, theoretical solutions are known and can 
be used. These ideas are illustrated in this section. 


The state vector Xj^ is transformed to a modal- amplitude vector c^^ 
real form of the eigenvector matrix as follows: 


with the 


X„ = 

m mm 


(3) 


Inserting equation (3) into equation (1) yields the modal equation of motion 


c 

m 


Ac + G A + Tu ~^D Ix^ 
mm me '-m m-'d 


(4) 


where 


A = U ~^A U 
m m m m 


Gm = 


u. 


- 1 , 


m 


in 


(5) 

( 6 ) 
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Matrix A is composed of the system eigenvalues in block diagonal form as 
m 


X 


m. 


1 


0 


A 


m 


0 


m,i 


• X 


m,n 


( 7 ) 


The block elements X^ ^ are assumed to be distinct eigenvalues, either real or 
complex-conjugate pairs. For real eigenvalues. 


X . = a. 
m,i 1 


and for complex eigenvalues. 


X . = 
m, 1 


a. 


X 


0 ). 

X 


where and o)^ are the real and imaginary parts. 

The discrete solution to equation ( 4 ) at time of 


+ T 

m o 

for period T is 

mm mmo mmco '-mdom-'do 

where matrices ^ , T , and G have the form 
mm m 


^ = 
m 




(j) . 

m, 1 


m,n 


r = 

m 


'm, 1 . 




ra,n 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


( 13 ) 
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and 


G = 
m 


1 






The solution and form for (b are shown later* 

iti 

The solution for ^ is 

m 


^ - exp (A T) 
m ^ m 


with block element 


K -! " exp(a.T) 
m, 1 1 


for X . real and distinct and 
m, 1 


(b , - exp(a.T) 
m, 1 


cos o).T 

X 


-sin o).T 

X 


sin o).T 

X 


cos 0 )^T| 


for X . complex. 
m,x 

The solution for F is 

m 


J T 

"" - t)] d-c 

. 


The block element is 


’m 


/I \ 

,i “ laJ 


(exp(a^T) - 1) 


for X . real and distinct and is 
m, X 


'm,x 


. 2 ^ 2! 

la. + 0 ). 

X X 


1m/i 

^12m,i 


^12m,i I 

1m, i I 


( 14 ) 


(15) 


( 16 ) 


(17) 


(18) 


(19) 


( 20 ) 
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( 21 ) 


for \ . complex / where 

m, 1 


, = (exp(a.T) ) (a. cos co.T + co. sin co.T) - a. 
M1m,x ^ X X XX X X 


. = (exp(a.T) ) (a. sin w.T - co. cos w.T) + o). (22) 

^12m,x XX XX X X 

The modal equations are now in explicit decoupled form and can now be used to 
design stability feedback gains- Putting equations (12), (13), and (14) into equa- 
tion (11), neglecting the disturbance term since it does not affect stability gain 
design, and partitioning the equations relating to the unstable flutter mode yields 
the block equation 


(23) 


The control-feedback signal is related to the modal amplitude c^^ ^ by the 

control-feedback gain matrix as 




Do Cm i (to) 
c m, X o 


(24) 


where 


°c = [’'Ic ’'2 c] 


(25) 


since there is one control surface and two modal amplitudes relating to the complex, 
unstable flutter mode. In this paper, a pole-placement technique is used; the deri- 
vation of the feedback gains for the pole-placement approach is given in 
appendix B- 

The closed-loop equation for the controlled mode is determined by substituting 
equation (24) into equation (23) as follows: 


c 


. (T ) 
m,x m 



m 


i^m, 


.D 
X c 


Ic .(t ) 
J m,x o 


(26) 


For stability, the eigenvalues of the 2x2 closed-loop matrix must have a magnitude 
less than unity. It is assumed that all except the term relating to the flut- 

ter mode, have open-loop eigenvalues located within the unit circle. This design 
procedure assumes that higher order modes are influenced as little as possible. 
Stability of the total closed-loop system is described in a later section. 
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ESTIMATOR MODEL FOR WIND GUST AND EXTERNAL- DISTURBANCE STATES 

If the modal amplitudes ^ of equation (23) could be measured directly/ the 

stability problem would be solved. Unfortunately/ only physical states can be mea- 
sured/ and these states are related to all (eq. (11)). Therefore/ an estimator 

is required to obtain the modal amplitudes in the feedback configuration described 
previously. 

During the period of this research/ it became obvious that the wind- gust veloc- 
ity Wg and the external- disturbance state vector must be estimated to get a 

good approximation of the two unstable modal- amplitude coordinates. One approach is 
to augment the states in equation (2) with and Wg and then solve the modal 

equations as described previously. The approach used in this paper is to augment the 
modal equation (eq. (11)) after solving for the discrete solutions for Wg and X^ 
and the effect of W^ and X^ on given as matrix (|ji^ in that equation. 

The wind gust is assumed to have a first-order/ slowly varying time response 
with time constant t defined as 

g 

w = - ~ W (27) 

g g 

g 

with the solution 


W (T ) = exp(-T/T )W (t ) (28) 

g m ^ g g o 

The differential equation representing the external disturbance is 


= Vd Vg 


(29) 


where A^ and are the external-disturbance state matrix and forcing-function 

input matrix. The solution/ for a step input W^ applied to a zero-order hold/ is 


^d<V = exp(A^T)X^(tQ) + Aa-‘'[exp(A^T) - l]B^Wg(to) 


(30) 


for A^ nonsingular. The two eigenvalues of A^j are defined as and 

are used in the solution for (1>^. For simplicity, we can define 


= exp(A^T) 


(31) 


and 


fa - 


(32) 
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and substitute into equation (30) as follows: 


*d'V = 




r.w (t ) 

d g o 


( 33 ) 


The general solution for (b in equation (11) results from solving the forcing- 
^ m 

function disturbance term in equation (4) as follows: 

T 

[ W’=o'’'n.> V ■ dT (34) 

o 

The solution to this equation requires the use of equations (A21), (A23), (28), 

and (33) and the solution to three convolution integrals. The derivation is given in 

appendix C. For the text of this paper, (|;^ is defined as 


(})=[(1j(1;] (35) 

^m g d 

where (\> and c|; relate the initial wind-gust velocity and external-disturbance 
state to^the modal-coordinate states. Using these augmented states from equa- 
tions (28), (33), and (35) and inserting them into equation (11) yields the equation 
for the discrete modal state vector X^: 


X (T ) = A X (t ) + 
o m o o o 


B A (t ) + L ^(t ) 
o c o o o 


where 





m 




A - 
o 


exp(-T/i:^) 


B = 
o 


r G 

m m 


0 

0 




0 



(36) 


(37) 


(38) 


(39) 
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and Lq is a vector relating the external wind ^ to the state vector X^. For 
this paper, is defined as the middle-column vector in equation (38) . The aug- 
mented states Wg and also cause changes in the form of equation ( 1) . By 

making use of equation (3), the acceleration at time T^^ is 


X (T ) = [a U W 
mm *- m m m 


F lx (T ) + B A (t ) 
m-* o m m c o 


(40) 


where and F^ are part of the derivation given in appendix C. 

The measured acceleration output , assuming negligible sensor dynamics, 

is defined as 


Y (T ) = C X (T ) + D A (t ) + 0 (T ) 
om oom oco om 


where 0 is measurement noise, 
o 


(41) 


Co=«o[%«m ^m] 


(42) 


^o = 


G B 
o m 


(43) 


and is a transfer matrix with the number of rows equal to the number of measure- 

ments. Each row of G^ has all 0*s and a 1 located in the column corresponding to 
the state that is being measured. 


KAIiMAN-FILTER MODEL 

The discretized model of the system dynamics and output measurement is described 
by equations (36) and (41) and is shown in the upper portion of figure 3. Initial 
conditions for this model are explained below. The initial state vector X^(0) is 
assumed to be a random vector with Gaussian distribution and has an expected value 


e{Xq(0)} = X^(0) = Xq(0) 

A 

where X^ is the best estimate of the initial state vector. 
The covariance P is defined as 

covFx (0); x'^(O)] = P(0) > [0] 
o o ■* 


(44) 


(45) 
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where the variance X is 

o 

X = X - X (0) 
o o o 

With similar notation for other variables, the Gaussian, discrete white plant 
noise ^ is 


K{?(t)} = 0 


and the covariance is 


cov[^(t); 


tT 


where ^ is the plant-noise intensity and 6. is the Kronecker delta defined 

tT 


6 . = 1 
tT 


(t = t) 


6 ^ = 0 
tT 


(t ^ t) 


Similarly, the Gaussian, discrete white measurement noise 0 is 

o 


E{e^(t)} = 0 


and covariance is 


covfe (t); 9 = ? 6 


O tT 


where 9 is the measurement-noise intensity. In addition, X (0), ^(t) , and 

o ^ 

0 ( t) are mutually independent for all t and t- 
o 

With equations (48) and (51), the state process noise R is defined as 


n T 

R = L 

o o 


and the variance of measurement noise is defined as 

Q. = 
o o 


Both R and are assumed to be constant during a run. 


( 46 ) 


(47) 


(48) 

as 


(49) 


(50) 


(51) 


(52) 


( 53 ) 
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with the above definitions and the definition for conditional probability/ the 
equations for the predict cycle are 


X (T t ) = A X (t t ) + B A^(t ) 
o m o o o o o o c o 


(54) 


P(T t ) = A P(t t )A + R 
mo o o o o 


(55) 


A A A 

where X^ represents the estimated value of X^ and A^ and represent approx- 
imate matrices for A^ and The approximated matrices have been used both for 

off-design cases and for cases with less than full-state estimation. 

The equations for the update cycle are: 


G, = P(T It )C *^rc P(T It )C + Q 
k m!o o'-o mlo o ^oJ 

(56) 

Y (T ) = C X (T It ) + D A (t ) 
om oomlo oco 

(57) 

X (T It ) = X (T It ) + G, (y (T ) - Y (T )] 
omim omio k^om om^ 

(58) 

P(T It ) = P(T It ) - G, C P(T It ) 
ml m ml o k o ml o 

(59) 


where is the Kalman gain niatrix. The structure for the discrete Kalman-f ilter 

model is shown in the lower portion of figure 3. 

The one remaining link is the relationship between the estimated state vector 
Xq and the control- feedback signal A^. The equation relating these variables is 


= 

c o 


D G X 
c c o 


o o 


(60) 


where G^ is a matrix with two rows of elements that are all 0 except for a 1 in 
each row located in the columns that allow the two modal- amplitude estimates of the 
unstable mode to be controlled. 


Closed-Loop System 

A closed- loop stability check is made by looking at the stability of the Kalman 
loop by itself and by looking at the complete closed-loop control system that 
includes both the plant and the Kalman model. ^By using equations (41)/ (54)/ (57)/ 
and (58)/ the updated estimated state vector X^ becomes 
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( 61 ) 


X„(T„It„) = Ta^ - GvC„A„]x„(t^|t„) 
o iti'ni Lo koo-io oio 

[»o - '=kt=o®o * “o - folK'^o’ 
* Vo>'o'''n,' * WV 


Stability of the Kalman loop requires that the eigenvalues of [A^ - ^k^o^o ' 
within the unit circle- The state equations for the complete discrete closed-loop 
system are 


where 


o m 


X (T T ) 
o m ' m ' 


'11 


21 


12 


22 


o o 




X (t t ) 
o o ' o ' 


G, C L 
k o o 


c(to) + 


0 (T ) (62) 

o m 


®11 " 


(63) 


12 O C C 


®21 “ %Po^o 


(64) 

(65) 


S~~ = A - G.CA + BDG - GiTcB - CB + D - dIdG 
22 o ^ o o o c c o o o o o o -i c c 


( 66 ) 


Stability requires that the eigenvalues of the matrix 
unit circle. 


11 12 

5- - 

21 22i 


lie within the 


Suboptimal Filter 

The optimal Kalman filter for a plant produces a time- varying gain in which the 
entire state is estimated and sufficient observations are processed to ensure that 
all the required states are observable. Moreover, in the presence of random, unknown 
error sources, process noise must be added to maintain an adequate estimate of the 
state uncertainty. Since this requires the propagation and correction of the state 
covariance matrix for each observation, it is not always practical to use an optimal 
filter in a real-time control system. The usual alternative is to simulate the opti- 
mal process off-line until the Kalman gain stabilizes, and then use the steady-state 
Kalman gains as a constant-gain complementary filter. However, even this approach 
requires the updating of a large state vector and the processing of a large number of 
observations. It thus becomes pertinent to seek a smaller state vector and fewer 
observations which could be used to generate a smaller dimension, constant yet stable 
gain. Such a system is defined to be suboptimal, and the equations used for the 
suboptimal system in real time are equations (54), (57), (58), and (60). The off- 
line calculations for the optimal Kalman gain are made assuming full-state esti- 
mation, with Aq, Bq, Cq, and Dq equal to A^, B^, C^, and D^. 
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The results for two different methods of reducing the number of estimated states 
are presented in this paper. The first method is to partition the matrix after solv- 
ing for the optimal gain Gy^* In particular, the method involves calculation of the 
reduced matrices by partitioning both rows and columns of A^, partitioning only rows 
of partitioning only columns of C^, and leaving unmodified. The second 

method involves calculation of the reduced Gy, directly. This is accomplished by 
using A^ and in equations (55) , (56) , and (59) . Results for both methods are 

described hereinafter. The actual states used are presented in the "Results” 
section. 


RESULTS 


This section presents the results for two different sets of closed- loop design 
eigenvalue (pole) locations for the flutter mode. For each set of pole locations, 
the statistical accuracy of several state variables resulting from wind-gust distur- 
bances is plotted as a function of dynamic pressure. Plots are presented for the 
case in which perfect modal feedback is used, for the case in which full-state esti- 
mation is used, and for a few cases in which reduced- state estimation is used. In 
addition, some typical data signals and step responses are presented. Stability of 
the filter and total closed-loop system has been checked with equations (61) and (62) 
for every case considered. 


The system defined by equation (1) has 33 states: has 11 states defined by 

5 bending states, 5 torsional states, and 1 control- surface state; X has 11 states 
defined by the derivatives of the X^ states; X^ has 10 states, since the 
unsteady-aerodynamics terras are modeled by a second-order differential equation for 
each wing panel; and X^ has 1 state. The unsteady aerodynamics is based on Jones* 
theory for two-dimensional incompressible flow (ref. 14) and is a good approximation 
to the lift-deficiency function. State X^ is always included, but the compensator 
can be effectively nullified by setting the numerator and denominator frequencies 
approximately equal and large compared with the flutter-mode frequency. Figure 4 
contains a plot of the open-loop eigenvalues as a function of dynamic pressure for 
the 10 flexible-wing modes. The first flexible mode is the flutter mode, with flut- 
ter onset occurring at approximately 54.0 kPa. All modes have a structural damping 
ratio of 0.03. Table I contains the open-loop eigenvalues at the design dynamic 
pressure q of 82.7 kPa. The eigenvalues are numbered for reference. The equiva- 
lent open-loop eigenvalues are also shown for the discrete plane with a sampling 
period of 0.005 sec, which is the period used for all runs in this paper. For this 
sampling period, modes 19 and 20 have a frequency greater than one-half the sampling 
frequency. 


Stability is the most important aspect in flutter control, and it must be main- 
tained throughout the flight regime. A second important aspect is the wing response 
and control-surface activity, resulting from a random wind-gust input disturbance 
(eqs. (A29) and (A30)). The approach used in this paper, for each of the variables 
recorded, is to normalize the rms value about a zero mean by the rms value about a 

zero mean of the wind gust c, • The recorded standard deviations are for the 

W 

g 


wing-tip bending 



wing-tip angular rotation 


, control-surface deflection 
“5 


Og, control- surface angular velocity 

a . 

u 


ot, and the 
o 


control-torque 


input command 
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TABLE I 


OPEN- LOOP EIGENVALUES 


[q = 82-7 kPa] 


Mode 


No 


Flexible 

wing 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Continuous-plane 

eigenvalue 


Real 

Imaginary 

15.98 

47.1 

15.98 

-47.1 

-38.41 

61.4 

-38.41 

-61.4 

-5.42 

139.8 

-5.42 

-139.8 

-6.31 

214.0 

-6.31 

-214.0 

-8.68 

263^5 

-8.68 

-263.5 

-8.93 

306^1 

-8.93 

-306.1 

-11.19 

379.8 

-11.19 

-379.8 

-13.93 I 

450^1 

-13.93 

-450.1 

-17.34 

605.8 

-17.34 

-605.8 

-21.10 

729.5 

-21 .10 

-729.5 


Discrete-plane 

eigenvalue 


Real 

Imaginary 

1 .0530 

0.2527 

1.0530 

-.2527 

.7867 

.2493 

*7867 

-.2493 

.7449 

.6264 

.7449 

-.6264 

.4653 

.8499 

.4653 

-.8499 

.2399 

.9270 

.2399 

-.9270 

.0384 

.9556 

.0384 

-.9556 

-.3050 

.8951 

-.3050 

-.8951 

-.5864 

.7253 

-.5864 

-.7253 

-.9111 

.1030 

-.9111 

-.1030 

-.7872 

-.4361 

-.7872 

.4361 


Actuator/ 

control 


surface 


21 

22 


-357.6 365.1 

-357.6 -365.1 


-0 


.04216 0.1619 

.04216 -.1619 


Unsteady 

aerodynamics 


23 

-72.11 

24 

-71.80 

25 

-70.27 

26 

-69.18 

27 

-33.16 

28 

-10.96 

29 

-10.93 

30 

-10.84 

31 

-10.71 

32 

-5.78 


0.6973 

.6987 

.7037 

.7076 

.8472 

.9467 

.9468 

.9473 

.9479 

.9715 


Wind 


gust 


33 


-0.6992 


0.9966 


External 

disturbance 


34 -31.1 

35 -240.0 


0.8559 

.3011 


Compensator 


36 


-500.0 


0.08208 

















Il 


Two sets of closed-loop pole locations are evaluated. One location is for a 
low-damped control mode = 0.14) and the other is for a higher damped control mode 
(C = 0.62); both pole locations have approximately the same damped natural frequency. 
The two design locations are given in table II for both the continuous and discrete 
planes. The relationship between the continuous-plane and discrete-plane eigenvalues 
is given in appendix B. 


TABLE II.- DESIGN POLE LOCATIONS FOR CONTROL MODE 


Set 

Continuous-plane location 

Discrete-plane location 

Damping ratio 

1 

-10 ± 70j 

0.89 ± 0.33j 

0.14 

2 

-58 ± 74j 

.70 ± 0.27j 

.62 


Four different cases have been simulated for each design pole location. The 
first case, which serves as a reference case for all other data runs, assumes full- 
state feedback to get a perfect estimate of the unstable modal- amplitude coeffi- 
cients. In addition, the control- feedback gains D^ are recalculated for each 
dynamic pressure q in order to maintain the controlled mode at the design pole 
location. Figures 5 and 6 show the feedback gains used to maintain the controlled 
mode at the design eigenvalue. For all other runs described in this paper, the 
feedback gains are kept constant with the values at q = 82.7 kPa, shown by the 
dotted lines in figures 5 and 6. 

The other three cases use three output measurements: wing- tip acceleration 

(grid point 5 in fig. 1), control-surface position, and control-surface angular 
velocity. All feedback numbers that would be manipulated in a computer (estimator 
model and stability feedback gains) are maintained constant over the complete range 
of q. 


The second case is for a full-state estimator design (X^ in eq. (37)). For 
this case, the open- and closed- loop eigenvalue variations of the first two modes are 
illustrated in figures 7 and 8. These figures are shown in the continuous plane for 
ease of understanding. For the design eigenvalue at 0.89 ± 0.33j, the mode 1 closed- 
loop eigenvalue shows a small variation, although damping is very small at low values 
of q. The mode 2 closed- loop eigenvalue illustrates a larger variation than the 
open-loop case for off-design values of q. With the design eigenvalue of 
0.70 ± 0.27j, both the mode 1 and mode 2 eigenvalues illustrate a wide variation over 
the range of q. 

The last two cases are for reduced-state estimators (suboptimal filter) and 
represent different states for the two design pole locations. The partitioning 
method is used whereby rows and/or columns of nonestimated states are removed from 
the full-state estimator design. The goal for reduced-state estimation is to reduce 
the number of computer calculations while maintaining a stable system over the range 
of q. Table III has the total number of estimated states and the eigenvalue number 
as given in table I for each design eigenvalue. 
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TABLE III.- REDUCED-STATE ESTIMATION 


Design eigenvalue 

Case 

Total no. of states 

Eigenvalue no. 

0.89 ± 0.33j 

3 

13 

1 to 8,21,22,27,33,34 

.89 ± 0.33j 

4 

9 

1 to 6,21,22,27 

.70 ± 0.27j 

3 

17 

1 to 12,21,22,27,33,34 

.70 ± 0.27j 

4 

15 

1 to 12,21,22,27 


Notice that fewer flexible-wing modes are required to maintain stability for the 
low-damped control mode. A possible reason is that the modal amplitude of the lower 
damped mode is more dominant compared with the higher order modes, making the lower 
mode more observable in the output measurement compared with the higher damped con- 
trol mode. Both simulation and analysis have shown that the actuator/control-surface 
mode must be estimated to maintain stability, while only 1 of the 10 unsteady- 
aerodynamics modes requires estimation. The mode selected has the largest Kalman 
gain feedback values of all the unsteady-aerodynamics modes. For stability it is 
possible to either use or ignore the wind-gust estimate and the external-disturbance 
states. 

Simulation results for the four cases described are illustrated in figures 9 
to 18. The first five figures are for the closed-loop design eigenvalue at 
0.89 ± 0.33j and the last five figures are for the design eigenvalue at 
0.70 ± 0.27j. Numbers on each curve correspond to the appropriate case- 

Two anomalies are readily observable. The first is that curve 2 shows less 
bending and less torsional sensitivity to wind-gust disturbance than does curve 1 at 
high dynamic pressures. One possible explanation is that because the estimator for 
curve 2 is not optimized for the off-design cases it is probable that the increased 
modal interaction results in modal contributions that subtract from the overall 
total. It is just as likely that these contributions add to the total, resulting 
in a larger state variation. The second anomaly is that at low values of q the 
control-surface activity is greater for the reference case than for the other three 
cases. The reason may be that the low-pass filter action in the estimator results in 
a less sensitive control system than the case with perfect modal feedback. 

As expected, cases 3 and 4 show greater activity for rms wind gust than the 
first two cases. Comparison of the cases shows that less rms variation is obtained 
by estimating the wind gust and one of the two external-disturbance states. The 
final data point at q = 96.5 kPa is not shovm for curve 4 in figures 9 to 13 since 
both the closed-loop eigenvalues and the simulation data showed an unstable system. 

In figures 14 to 18 the closed- loop system was barely stable at q = 96.5 kPa, and 
the simulation showed a large increase in the accuracy of each variable. 

Several simulation runs were made for cases with measurement noise. A random 
noise with a standard deviation of 10 percent of the measurement was used for the 
accelerometer and control-surface angular velocity and a noise with a standard devia- 
tion of 1 percent was used for the control-surface position. In addition, a bias 
error of 0.025g was used for the accelerometer and a bias error of 0.005 rad/sec was 
used for the control-surface velocity. The simulation runs had statistical data very 
close to that shown for the perfect-measurement cases. 
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Some typical simulation plots are shown in figure 19 for the case of full-state 
estimation- The dynamic pressure was 82-7 kPa and the design eigenvalue was 
0.89 ± 0.33j with random and bias measurement errors. The plots in the first part of 
figure 19 represent the measurements of wing-tip acceleration h^ , control-surface 
angular velocity 6 , and control- surface deflection 6, and wind-gust velocity W^. 
The second part of figure 19 contains the wing- tip bending h^, wing-tip angular 
rotation wing-tip velocity h^, and wind-gust velocity W^. The top two curves 

in the third part of figure 19 are the actual modal-amplitude coefficients of the 
controlled mode while the bottom two curves are the estimated modal amplitudes that 
are used for control. 

Figure 20 shows some typical step responses resulting from an external command 
applied at the same input as feedback signal for q = 82.7 kPa with a design 

eigenvalue of 0.89 ± 0.33j. The four cases illustrated are perfect modal estimation 
(fig. 20(a)), full-state estimation (fig- 20(b)), 13-state estimation (fig. 20(c)), 
and 9-state estimation (fig. 20(d)). There are seven plots in each figure; h^, 
a^, 6, two plots of estimated modal amplitude, and two plots of actual controlled 

modal amplitude. In figure 20(a), the estimated modal amplitudes are set identical 
to the controlled modal amplitudes to get the perfect feedback. For full-state esti- 
mation, the estimated modal amplitudes are a very good representation of the con- 
trolled modal amplitudes, but as the number of estimated states decreases the quality 
of the estimated modal amplitudes also decreases. The noisy estimates cause 
increased noise on the control-surface position. The basic responses of the various 
states are approximately the same for all the cases, although there is a slight 
attenuation and a little oscillation in the h^ curve for the reduced- state estima- 
tion cases. 

A final set of data was taken for another suboptimal-f ilter case in which the 
Kalman gain is calculated directly using the reduced matrices A^ and C^. The 
results from this approach are generally not as good as the results from the parti- 
tioning method. For the reduced- state estimation cases for a design eigenvalue of 
0.89 ± 0.33j, all data are slightly worse than the data shown in figures 9 to 13. 

The closed-loop system is unstable at q = 96.5 kPa for the 13-state estimation, and 
the system goes unstable at q = 89.6 kPa for the 9-state estimation. For the 
reduced-state-estimation cases for a design eigenvalue of 0.70 ± 0.27j, the data 
indicate slightly better results for the 15-state-estimation case but worse results 
for the 17- state-estimation case. 


CONCLUSIONS 

A modal-control design approach has been used to stabilize a single-control- 
surface, unswept-wing model subject to bending-torsion flutter. The wing has five 
panels with five bending and five rotational degrees of freedom in addition to the 
one control- surface degree of freedom. The complete model also includes unsteady 
aerodynamics, wind gust, and external- disturbance states relating the wind gust to 
forces and moments on the wing. 

The modal analysis shows the relationship between the physical and the modal 
domains and the form of the matrices in the modal domain. Since the modal equations 
of motion are decoupled with block elements of one or two states, analytical solu- 
tions can easily be used to obtain feedback gains. The modal approach used in this 
paper is to partition the two corap lex- conjugate equations representing the flutter 
mode and the design feedback gains to stabilize this mode. Only two numbers 
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are needed to create a stable system. A procedure for calculating these gains for a 
discrete system for the pole-placement approach used herein is described in 
appendix B. 

The equations for a fixed-gain discrete Kalman filter with a reduced number of 
output measurements are described. The analysis includes equations Eor the estimator 
stability and for the stability of the total closed-loop system. Two approaches for 
using a suboptimal filter are included. 

Statistical data runs using a random wind-gust disturbance have been made for 
two different closed-loop design-eigenvalue locations for the controlled mode. One 
eigenvalue location is for a low-damped case (0.14 damping ratio) and the other is 
for a higher damped case (0.62 damping ratio). With three output measurements (wing- 
tip acceleration, control-surface position, and control-surface angular velocity) , 
the full-state estimator (36 states) stabilizes the wing over a dynamic-pressure 
range that is 80 percent greater than that of the uncontrolled system. The design 
point is at a 50-percent increase in dynamic pressure, and the feedback gains are 
maintained constant for all off-design dynamic pressures. 

Stability can be maintained with fewer estimated states for the low-damped 
eigenvalue case than for the higher damped case. For the low-damped eigenvalue case, 
a 13-state estimator model allows stability to be maintained over the entire dynamic- 
pressure range with only a small decrease in performance, while a 9-state estimator 
model allows stability to be maintained over most of the range (68-percent increase 
in dynamic pressure). Simulation runs were made for both the 17-state and the 
15-state estimator models for the higher damped eigenvalue case. Both models allow 
stability to be maintained over the complete dynamic pressure range with a small 
decrease in performance. The best performance results are obtained in the cases 
where wind gust and an external-disturbance state are estimated. Random measurement 
errors with a standard deviation of 10 percent of the measurement for both wing-tip 
acceleration and control-surface angular velocity and 1-percent standard deviation 
for control-surface deflection, along with small bias errors, appear to have negligi- 
ble effect on all of the runs described. 

All the data shown for the reduced-state estimators are for the approach where 
the reduced states are obtained by partitioning rows and/or columns of the matrices 
in the full-state estimator design. The partitioning method shows better stability 
and generally less state variations due to wind gust than the approach where the 
reduced Kalman gains matrix is calculated directly using a reduced plant model. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 26, 1982 
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APPENDIX A 


MATRICES FOR DEVELOPMENT OF SUBSYSTEM MODELS 

This appendix contains the definitions and forms of all matrices used to develop 
the structural model/ the aerodynamics model, the external-disturbance model, the 
wind-gust model, and the compensator model. In addition, a derivation for the 
structural-damping model is included. Most of the data for the structural matrices 
and the steady-state aerodynamics influence coefficients come from reference 13. 


Structural Equation of Motion 

The state vector for the wing is defined to have 11 states: 


X = 
w 


w 


w 


(AD 


where h^ is the 5-state bending displacement with positive axis down, a is the 
5 -state torsional rotation with positive leading edge up, and 6 is the trai ling- 
edge control- surf ace displacement corresponding to trailing edge down. (A list of 
symbols used in the appendixes appears after the references.) The equation of motion 
for the wing structure is 


MX+CX+KX=Gu+L+L^ (A2) 

sw sw sw s a d 

where M^, Cg, and represent the structural matrices for the mass-inertia 

terms, damping terms, and wing bending/torsional-stif f ness terms, represents the 

control-torque input matrix, u represents the control-torque input variable, L^ 
represents the aerodynamic lift and moments, and represents the external- 

disturbance forces and moments. The individual matrices are defined below. A devia- 
tion of Cg is included at the end of this appendix. 

The structural mass-inertia matrix Mg is defined as 



where 


= 




®65^65 


(A3) 


(A4) 
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«5 = 


” 65''^65 


(A5) 


and is a diagonal mass matrix, is a diagonal inertia matrix, is a 

diagonal mass-unbalance matrix, Sg is a control-surface mass-unbalance vector, Hg 
is a control-surface moment of inertia vector, is the moment of inertia about 

the control hinge, and V« defines the wing panel containing the control surface as 

o 


Vj - [0 0 



(A6) 


The wing structural-stiffness matrix is defined as 


K = 
s 




0 

0 


'^ 6“6 


(A7) 


where is the bending-stiffness matrix, is the torsional-stif fness matrix, 

and cOg is the control-surface structural radian frequency. Both and are 

full matrices and are calculated by taking the inverse of their respective flexibil- 
ity matrices. 


The control-torque input matrix G„ is defined as 

s 

fol 



1 


(A8) 


The matrix has one column relating to the single control surface being used. In 

this analysis, the control is modeled as a torque. The actuator torque is applied to 
the control-surface rotation axis, and an equal and opposite torque is applied 
directly to the corresponding wing section. 


Aerodynamics Model 

The unsteady-aerodynamics model contains Jones' exponential approximation of 
Wagner's indicial loading function for incompressible two-dimensional flow (refs. 14 
and 15) as a scalar multiplier of the three-dimensional steady-state lift and moment 


. 


20 



APPENDIX A 


distribution for the wing. Jones assumed a time-response function as 

follows: 


Ka(t) = 1 - 0.165 e^ 


b ) 


(A9) 


where V is the airspeed and b is the semichord. The unsteady-aerodynamics 
states are functions of both the wing-position states X and wing-velocity 

3L. W 

states as follows: 

w 


X = AX +B, X +B^X 
a aa law 2a w 


(A10) 


In return, X^ causes changes in as follows: 

cl cL 


L = CX +D^X +D^X 
a a a la w 2a w 


(All) 


where x^ is defined to have two five-element states and ^2a’ 


X = 
a 


la 


2a 


(A12) 


since K^(t) is shown to have a second-order response, 
a 

The aerodynamics matrices are defined as the following: 




AL 


D = -ff-) 

EA 

EAL 

k.EALV. 

la \2V/ 

-T 

-T 

6 6 

, -T 


C A 

C AL 

k.C ALV- 

0 0 


(A13) 


- -(f 


EA 




k_EAV. 
o o 


-T -T 

C A kX AV„ 

0 0 


(A14) 
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C = (q) 

ct 


EA 


T 

C A 


7^ V\^ /^0.1080S\ 

^0.006825 -j I ( b j 






0 0 


I '• ’'6“6 


= 2 . - <''> 


0 0 


0 I 


= 




-0.01365 ^ I 

I bV 


(-0.3455 ^II 


where 


-T T" 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 



2 


(A20) 


and where A is the influence-coefficient matrix for the steady-state aerodynamic 
lift distribution/ E is a diagonal matrix with elements equal to the distance 
between the elastic axis and the aerodynamic center at the ith station, is a 

diagonal matrix with elements equal to the distance from the center of pressure of 
the control surface to the hinge, L is a diagonal matrix with elements equal to the 
distance from the three-quarters of the chord to the elastic axis, k. is the ratio 
of control- surface lift slope to wing-section lift slope at the wing section contain- 
ing the control surface, q is the dynamic pressure, p is the density, I is an 
identity matrix, V is the airspeed, and b is the wing semichord. 


External-Disturbance Model 


The time response for the external disturbance has a form similar to that 

for the unsteady aerodynamics. The model used is the Kussner function (ref, 16) for 
the indicial response to a step of vertical gust velocity defined as 


K^(t) 



(A21) 
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In a manner similar to the wing unsteady lift and moment, the two-dimensional 
response is used as a scalar multiplier to the steady-state three-dimensional 

aerodynamic lift and moment distribution. As shown in equation (29), the differen- 
tial equation for the external-disturbance state vector is a function of the 

wind-gust velocity W^.: 


X, = A,X^ + B^W 
d d d d g 


with 




(A22) 


Since is a second-order response, X^ has two states: 




Id 


X 


2d 


(A23) 


The matrices for external-disturbance states are the following: 


Cd = (q) 


A 

EA 


-T 
C A 


rf) 

1 

1 

1 

1 


0.13 ^ 
_ 


0.565 


(A24) 


^d = 


-0.13 


-1.13 I 


(A25) 


B 


d 


0 

1 


(A26) 


Wind-Gust Model 


A Dry den wind model is used to simulate W 


Dryden spectra 




for turbulence velocity (ref 


for the transverse direction. 
17) is 


The 


•w <“g> 

g 


2 -»■ ] 

lt[l + 


(A27) 
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where Qg is the spatial frequency, is the wind- turbulence scale length (533 m 

(1750 ft) in ref. 17), and Og is the rms wind-gust intensity. By using the rela- 
tionship between spatial frequency and temporal frequency w as 

g 


0) = Q V (A28) 

g g 


and by assuming a discrete random^ nxmber generator with sampling period and 

band-pass frequency 2n;/T , the equations of motion and matrices become 

g 


X = A X + B W (A29) 

g g g g c 


and 


where X_ has two states and W_ is the white noise of unit intensity. 

g 

The matrices for the wind-gust model are the following: 


(A30) 


0 


A = 

g 



1 



0l 


B = 

g 


c = 
g 



where 


Y 


g 


V 

L 

g 


(A31) 


(A32) 


(A33) 


(A34) 


Compensator Model 

There are occasions when in-line compensation, such as a lag circuit, a lead 
circuit, or a combination, is to be installed. A compensator has been installed at 
the input to the actuator and is driven by a control- feedback signal A^ as follows: 


X^ = A^X^ + B^A 
f f f f c 


(A35) 
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and 


U = CfXf + 


(A36) 


where A^, C^, and are compensator matrices and is the compensator 

state. Given a transfer function with a lead frequency and a lag frequency 

and a gain E^z the compensator terms become 




(A37) 


Cf = - o^2f 




(A38) 

(A39) 


Complete Model 

The equations described previously are combined in block diagram form represent- 
ing the complete open-loop system in figure 21. The key equations are given in the 
blocks and the equation numbers are given above the blocks. The measured output 
Yq is shown as a function of all the states (eq. (41)). Equation (A2) requires some 
rearranging to get it into proper form for analysis. Combining equations (All), 
(A22)z and (A36) with equation (A2) yields 


X =DX +AX +GX + G^X^ + G^X^ + B A 
w ww ww aa ff dd cc 


where 


D = M ~V-C + D. 1 

w s s la-" 


A = M ~V-K + D- 1 
w s s 2a-* 


G = M ~^C 
a s a 


M “’c. 

d s d 


G- = M "^G C_ 
f s s f 


B = M ~^G E- 
c s s f 


(A40) 


(A41) 
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With these definitions/ the matrices for equation (1) are 


A = 
m 


D 

w 

I 

0 


B 


w 

0 

2a 

0 


(A42) 


B = 
m 


(A43) 


D == 
m 




(A44) 


Derivation for Structural-Damping Matrix 
The homogeneous differential equation for the wing structure is 


MX + CX + KX = 0 
s w s w s w 


(A45) 


By using the modal transformation. 




(A46) 


The equation of motion in modal form becomes 


Z + DZ + AC = 0 
s s s s s 


(A47) 


where 


U 

s s s s s 


(A48) 
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A = U ~V “^K U (A49) 

s s s s s 

Equation (A47) is completely decoupled with each differential state equation of the 
form 


• 2 

Z .+2C .CO ,Z . + 0) .z . = 0 
S/X S/X s,x s,x s,x s,x 


(A50) 


where cOg ^ is the radian frequency for mode i and Cg,i structural damp- 

ing ratio for mode i. The matrix Ag is calculated by letting the damping ratio be 
0 in equation (A45) and solving for the eigenvalues. The matrix Dg is determined 
next by inputting approximate values for C • With equation (A48), Cg is calcu- 
lated as 
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DERIVATION OF FEEDBACK GAINS FOR POLE- PLACEMENT APPROACH 

This appendix presents the derivation and final equations for calculating the 
stability feedback gains for a discrete control system, in which the desired pole 
locations are either in the unit circle for discrete pole placement or in the stan- 
dard S-plane for continuous pole placement- The derivation is for a complex- 
conjugate plant with one control input A^. 

The modal equation of motion is given in equation (4): 


c (t) = 
m 


A c (t) 
m m 


G A (t ) 
me o 


and the discrete form solution for the unstable flutter mode is presented in 
equation (23) : 


c (T ) = (t> c (t ) + [y g ]a (t ) 
mm ^m m o 'm m-* c o 


Equations (24) and (25) define the control feedback as 


c o 


D C (t ) 
c m o 


where 


D 

c 




Finally, the closed- loop expression in equation (26) is 


c (T ) = [(() + 

mm 


m 


Y g 0 Jc (t ) 
'm m C-* m o 


In the preceding equations, subscript i is deleted for convenience. 

Subscripts m and c and the arguments T^ and t^ will be dropped in the 
remainder of this paper for simplicity. If the desired poles are located in the 

continuous plane, the first step is to calculate the equivalent poles Z in the 
discrete plane. Defining and Z in terms of the real and imaginary parts gives 

Z„ = exp(a.T) cos (w,T) (B1) 

R a a 


Z_ = exp(a.T) sin (w.T) (B2) 

la a 
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where and are the real and imaginary parts of and Zj are the 

real and imaginary parts of Z, and T is the sampling period. 


By defining y, and g in terms of their components, equation (23) becomes 




- 


'*’11 ‘*’ 12 : 


*11 




v!si 


'11 '12 


-Y 


12 '11 


(B3) 


With the definitions 


V 


1 


Yll^l ^12^2 


and 


(B4) 


= -^12^1 


^11^2 


(B5) 


the closed-loop expression in equation (26) is 


s 






*11 *12 


*'21 '''221 




(B6) 


where 


*11 = *11 ^iNc 


(B7) 


^2 = *12 ^ ^l’^2c 


(B8) 


^21 = -*12 Vic 


(B9) 


**22 *11 ^ 2 ^ 2c 


(BIO) 
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The closed-loop eigenvalues of equation (B6) are found by using 


Z - 


^11 


-4) 


12 


= 0 


(BID 
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Z - 


‘'22 


with the solution 


Z = 


2 ^ 2 ^^ 12^21 


(B12) 


Substituting equation (B7) into equation (BIO) yields 


26.. + v.k. + v_k. 

^11 1 1c 2 2c . 

2 = j. 


4 Ufv.k. + v_k- - 46..r6._ + v.k. - v.k. 1 
2 V' 1 1c 2 2c'' ^12'^12 1 1c 2 2c^ 


(B13) 


For complex-conjugate eigenvalues. 


Z = Zr ± z^j 


(B14) 


Equating the real parts of equations (B13) and (B14) yields 


ko = ;r-{2Z„ - 24)^ - ) 

2c R 11 1 1c 


(B15) 


For comp lex- con jugate eigenvalues/ the expression under the radical in equation (B13) 
is negative* Equating the imaginary parts of equations (B13) and (B14) and substi- 
tuting equation (B15) yields 


^1 

1c 


+ z/] - - Si) 




(B16) 


6*« + 6..-^v^ 

^12 v ^ ^12 2 
2 


For the two real eigenvalues, either distinct or equal, Z is expressed as 


^R1 


^11 '*’22 1 


- 2V^^1 -'^22) ^"^2^: 
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(B17) 
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or 


R2 


*11 ^22 1 fZ - a 2 

_ _ + _ - 4,^^) + 44 . 


12*21 


Adding equations (B17) and (B18) and using equation (B13) yields 


- 24)^^ - ) 

2c R1 R2 ^11 1 1c ^ 


Subtracting equation (B19) from equation (B18) and substituting for k 


2c 


"^R1^R2 


k = 
1c 


- *11^ ^ *12^ ♦ ♦itfSl * ^ ^R 2 - 




''12 v„ * 12^2 


(B18) 

(B19) 

yields 

(B20) 
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DERIVATION FOR TRANSFORMATION MATRIX RELATING EXTERNAL- DISTURBANCE 

STATES TO MODAL COORDINATES 

The Kussner function is given in equation (A21) of appendix A and the model 
given in equations (29) and (A23) . The solution to the differential equation is 
shown in equations (30) to (33). The differential equation (eq. (29)) 




+ B^W 
d g 


and the general solution (eq. (33) ) 


X^(T ) = 4)^X^(t ) + r W (t ) 
dm ddo dgo 


are repeated here with the following definitions: 


A = 


'**11d ‘**12d 


‘**21d ‘**22d 


r. = 


Id 


'2d 


(CD 


(C2) 


From equation 


(A24) and with the following definitions. 



(C3) 


k 


1d 



(C4) 


k 


2d 


-1.13 


V 

b 


(C5) 
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the external-disturbance force and moment vector is written as 




- - (^) 




Substituting equation (33) and combining terms gives 


^d'’’.’ ■ (v)(®dK*,d>'ld'’^o> * *2d>'2d<'o’ * 1'd"g"o'l 


where 


hd= "Nd^id" 


»>2d = -’'ld'*>12d - °-^’'2d'^22d 


Yd= -^Id^ld - °'^’'2dY2d 


With actual values the following are obtained; 


4>^d “ 0*065 ^(exp(\^^T) + ej^CX^^T)) 


= (0.065 exp(X,^T) + 0.5 exp(X^^T)) 


Id 


2d 


= 1 - 0.5(exp(X^^T) + expCX^^T)) 


The general solution for (J; is given in equation (34) as 


/.T 

C ^ -1 

(\> ,T )xJt 1=1 [4) (t - t)1[u D ]{x^(t)} dT 


Using equations (A44) , (A22), and the solution for in equation (A41) we 


r^m [m "^1 


dx 


(C6) 

(C7) 

(C8) 

(C9) 

(C10) 

(C11) 

(C12) 

(C13) 


get 

(C14) 
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Substituting equation (C7) yields 

-T 


in a 


r ^ 

,t ,T )X,(t ) = J [6 (t - •c)1{h }f(t>,.(i:)X (t ) 
o mdo ni la Ido 




where 


(M - {?)[”„■’ 


Rewriting equation (C15) gives 

4»„(X.,t ,T )X (t ) = 6 W (t ) + <|>. -X (t ) + (P~^X (t ) 
mdomdo ^ 9^0 Id Id o 2d 2d o 


where the three convolution integrals are 




J - 1^)) 


dT 


L^O 


{h } 


^1d 




(h } 


^2d 


1 


T 

m 


[$^(T)](4>2a(t - X)) dx 


{h } 


By letting 

^d = f‘*'ld ‘*'2d^ 

and with equation (35) 

(|> = [ <|> 4» a ] 

the discrete solution for the modal coefficients in equation (11) becomes 


c(T)=$c(t)+rGA(t)+(|>W(t) + 4»^X^(t ) 

mm mmo ramco ggo ddo 


(C15) 


(C16) 


(C17) 


(CIS) 


(C19) 


(C20) 


(C21) 


(C22) 
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The solution for each convolution integral 4^1 d' ^2d shown below* A 

typical block matrix for is given in equations (16) and (17) with elements 

defined in equations (Bl) and (B3). The equations for ”^d 

defined in equations (Cll) to (C13) . Combining all of these equations, the solution 
for equations (C18), (C19), and (C20) can be determined. The result below does not 
include the final multiplication by vector h^^, which would only affect the final 
gains. Since this is a typical block, subscript i is dropped for convenience. 

With the definitions 


= exp(X^^T) 

(C23) 

= exp(X2^T) 

(C24) 


the elements of each matrix are the following, with equations (C25) and (C26) relat- 
ing to equations (C27) and (C28) relating to equations (C29) and 

(C30) relating to (\) : 


‘*'11, Id ‘*’22, Id 


= 0.065 


V 



+ {cc- - e^) 

(a - -H 0," 

<^»12 ~ ^2d^^‘*’l1 ~ ^ 2^^ 

fa - 


(C25) 


4*. 


12, Id “‘*'21, Id 


„ i»i2^°‘ ~ Nd^ ~ ‘'^^'*’11 ~ ^ 1 ^ 

0.065 \ 99 

^ 1. ~ ^1d^ 

^ ‘*'12<^°‘ “ ^2d^ ~ ‘'^^‘**11 ~ ^ 2 ^ 
fa - . 0,^ 


(C26) 




11, 2d ‘*'22, 2d 


V /o-065[o) 4>,3 ~ Sd^^‘*'l1 ~ ^1^^ 

V (a - X^^) + 0) 

0.5[o 3(}>^2 + (« - >^2d^^‘*'l1 ~ ^2 ^^^ 


4>. 


= - 4 >, 


^ Jo.065[(|)^^(a - X^^) - a)f<t»^, - e^)] 


12, 2d ^21, 2d b 


fa - X^^f + 


0 . 5 [ 4 )^ 2 f “ “ ^ 2 d ^ " “(‘*’11 “ ® 2^^1 


(a - X^^l + 0) 


(C27) 


(C28) 
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'^22, g 


Y11 


0.5 



+ f« - 

(a - + (/ 




‘*'12, g ■‘^’21,g ^12 


+ 


0 5 ' Sal - m(»,, - 

■ I 

~ ^2d) ~ “Kl ~ °2il 

(« - J 


(C29) 


(C30) 


where and y^^ are defined in equations (21) and (22). 

Portions of the derivation presented here are also used to derive W and p , 
used in equation (40). Equation (1) is given as mm 


X = A X + B A + DX, 
m mm me m d 


Using equations (A23) , (A44), and the solution for in equation (A41) the distur- 

bance term becomes 


D 

m d 



(C31) 


Substituting equation (C3) and using the definitions in equations (C8) to (C10) 
results in 



(V*g * " 



We can now define 


(C32) 


(?)(V--)[*1d ♦2d] 


(C33) 


(C34) 
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SYMBOLS 


Symbols Applicable to Entire Paper 
control-feedback signal 
external-disturbance state matrix 
open- loop-model state matrix 

external-disturbance forcing- function input matrix 
open- loop-model control input matrix 
augmented discrete modal-control input matrix 
augmented discrete output modal state matrix 
modal- amplitude vector 
control-feedback gain matrix 
open- loop-model disturbance input matrix 
augmented discrete output- control input matrix 
discrete matrix relating external disturbances to 

transformation matrix relating estimated modal- amplitude vector to modal 
amplitude used for control 

Kalman gain matrix 

control input matrix for modal equation 

transfer matrix for output measurements 

ith element of matrix 

wing displacement for ith panel 

wing displacement for fifth panel 

identity matrix 

imaginary number 

control- feedback gains 

vector relating external random wind to state vector 
covariance matrix 


variance of measurement noise 



q 


dynamic pressure 


R state process noise 

1 '’"^22 submatrices for discrete closed-loop state matrix 

s Laplace transform 

T sampling period for control 

Tj^ time at end of sampling period 

t time 

t^ initial time 

eigenvector matrix 
u control-torque input 

Wg wind-gust velocity 

VI^ discrete matrix relating wind gust to 

unsteady-aerodynamics state vector 
a 

X^ external-disturbance state vector 

X^ compensator state vector 

X^ state vector for complete open- loop model 

Xq augmented discrete modal state vector 

X^ wing state vector 

measured acceleration-output vector 

real part of eigenvalue ^ 

a . wing torsional rotation for ith panel 

w,i ^ ^ 

wing torsional rotation for fifth panel 

vector defining discrete transformation between initial value of wind gust 
and external-disturbance state 

r control input matrix in the discrete modal solution 

Y . ith element of T 

*m,x m 

6 trailing-edge control- surface displacement 

6. Krone cker delta 

tT 
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m 


m, 1 


^1d'^2d 

m 

a 


T 

T 

9 


dcimping ratio 

Gaussian, discrete white measurement noise 

measurement-noise intensity 

open-loop eigenvalue matrix 

ith element of A 
m 

eigenvalues of 

Gaussian, discrete white plant noise 
process-noise intensity 
standard deviation 
time variable 

first-order time constant for wind-gust approximation 
discrete modal state transition matrix 


m 


m,i 


'I'd 


'I', 


m 


0) 


discrete transition matrix for external disturbance 


ith element of $ 


m 


submatrix of cb 

] 

submatrix of <b 


m 


m 


discrete transition matrix relating external-disturbance states to the 
modal-amplitude vector 


imaginary part of eigenvalue 

a. 

Superscripts : 

— mean value 

estimated value 
variance 

-1 inverse 

T transpose 

Special operators: 

e{ } expected value of a vector 


m 


exp 


exponential 


[ ] matrix 

cov [ ] covariance 

(x|y) conditional density of X given the value of Y 


Be 

Bf 

Bla'B2a 

b 


C 

C, 


s 

Bla'B2a 

E 


Symbols Applicable to Appendixes Only 

influence-coefficient matrix for the steady-state aerodynamic lift 
distribution 

unsteady-aerodynamics state matrix 
compensator state matrix 
wind-gust state matrix 
submatrix of A^^ 
submatrix of 

compensator-control input matrix 
wind-gust forcing-function input matrix 

unsteady-aerodynamics forcing-function input matrices 
wing semichord 

unsteady-aerodynamics output state matrix 
external-disturbance output matrix 
compensator-output state matrix 
wind- gust output matrix 
structural-damping matrix 
defined in equation (A14) 

diagonal matrix with elements equal to the distance from the center of 
pressure of the control surface to the hinge 

modal structural-damping matrix 

submatrix of 

unsteady-aerodynamics forcing-function output matrices 

diagonal matrix with elements equal to the distance between elastic axis 
and aerodynamic center at the ith station 

defined by equation (C3) 
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compensator-output-control input matrix 
submatrix of 
submatrix of 

m 

submatrix of 

m 

control- torque input matrix 
defined by equation (A5) 

control-surface moment of inertia vector 
defined by equation (C16) 
diagonal inertia matrix 

moment of inertia about the control hinge 
time response for Jones' approximation 
time response for Kussner function 
bending-stiffness matrix 
wing structural- stiffness matrix 
torsional-stiffness matrix 

ratio of control-surface lift slope to wing-section lift slope at the wing 
section containing the control surface 

constants 

diagonal matrix with elements equal to the distance from three-quarters of 
the chord to the elastic axis 

aerodynamic lift and moment vector 

external-disturbance force and moment vector 

wind-turbulence scale length 

matrix containing structural mass and inertia terms 

diagonal mass matrix 

diagonal mass-unbalance matrix 

defined by equation (A4) 

control-surf ace mass -unbalance vector 

sampling period for wind-gust random- niimber generator 



U modal transformation for free vibration 

s 

V airspeed 

Ve wing panel containing control surface 

o 

defined respectively in equations (B4) and (b 5) 
white noise of unit intensity 
Xg wind-gust state vector 

^1a'^2a components of unsteady-aerodynamics state vector 

Xid#X 2 d components of external-disturbance state vector X^ 

Zj discrete-plane imaginary root 

Zj^ discrete-plane real root 

Zg structural modal coordinate 

a, real part of desired pole 

d 

defined in equation (C10) 
band-pass frequency for wind model 
Yid'Y2d components of 

YiifYi 2 elements of matrix 

diagonal eigenvalue matrix for structure 
desired pole 
wind-gust intensity 

^11 '*^12 elements of the open-loop discrete modal transition matrix 

^11d''*’l2d''*'21d''{’22d components of 

^11 '*^12 '^21 '^22 elements of the closed-loop discrete modal transition matrix 

(})^ Dryden spectra 

g 

spatial frequency for wind-gust model 
6)^ imaginary part of desired pole 

CO temporal frequency for wind-gust model 
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conpensator frequencies 




0)g control-surface structural radian frequency 

One dot over a symbol indicates first derivative with respect to time. Two dots 
over a symbol indicate second derivative with respect to time. 
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Figure 1.- Simplified model of unswept aircraft wing 
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Figure 4.- Flexible-wing modes as a function of dynamic pressure 
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Figure 7.- Eigenvalue variations of first two modes as function of dynamic 
pressure for case of full-state estimation with mode 1 closed-loop design 
eigenvalue at 0,89 ± 0.33j (s = -10 ± 70 j) for q = 82.7 kPa. 
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Figure 16*- Root-mean- square control- surface rotation for 1 m/sec rms wind gust 
as a function of dynamic pressure for design eigenvalue of 0.70 ± 0.27j 
(s = -58 ± 74j) . 
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(a) Perfect modal estimation. 


Figure 20.- Step responses resulting from external command at input of A^. 
q = 82.7 kPa; Design eigenvalue = 0.89 ± 0.33j. 
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